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Abstract. —Nucleotide sequences of the chloroplast trnL-trnF intergenic spacer were obtained for 
21 of the 22 indigenous Asplenium taxa presently recognized from New Zealand. Nucleotide 
sequences of the chloroplast rbcL gene were also obtained from eleven New Zealand species 
representative ot the diversity found in the trnL-trnF intergenic spacer. Phylogenetic analyses of 
these chloroplast sequence data indicate that the Asplenium species of New Zealand are not 
monophvletic. More specifically, the Asplenium species participating in hybridization in New 
Zealand form a closely related ‘Austral’ group, whereas the non-hybridizing species have closer 
affinities to species from outside New Zealand. Within the Austral group, three well-supported 
sub-groups are recognized, represented by the species A. bulbiferum, A. flaccidum, and A. 
obtusatum. Dating analyses reject an 80 million year old vicariant origin for any of the Asplenium 
lineages in New Zealand, and the distributions of the many Asplenium species disjunct between 
New Zealand and elsewhere appear best explained bv long-distance dispersal. The likely 
chloroplast/maternal parent for each of the New Zealand octoploid species is discussed. 


Asplenium, with approximately 700 species worldwide, appears to have a 
complex history of hybridization and auto- and allopolyploidy, and is a model 
group for the study of fern evolution (e.g., Manton, 1950; Wagner, 1954; Lovis. 
1977; Reichstein, 1981). Recent DNA sequencing studies (e.g., Murakami et al., 
1999a; Yatebe et al., 2001; Gastony and Johnson, 2001; Pinter et al., 2002; Van 
den heede et al., 2002, 2003) have shed light on the evolutionary history of this 
iconic group of ferns. Here we provide chloroplast DNA sequences for 21 of the 
22 indigenous Asplenium taxa in New Zealand (Table 1; Brownsey and Smith- 
Dodsworth, 2000; Brownsey, 2003), where the genus is often ecologically con¬ 
spicuous and comprises about 10% of the indigenous fern flora. 

Early taxonomic studies of New Zealand Asplenium (e.g., Hooker, 1855; 
Martin. 1920) indicated hybridization was rife and found the different ‘kinds’ to 
grade together. In his account of New Zealand Asplenium, Allan (1961. p. 75) 
remarked that the species were “very ill-defined", and that while many species 
appeared to respond “markedly to environmental conditions ... [,] there is also 
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Tahi.k I. Indigenous Now Zealand Asplenium, together with additional species included in study. Ploidv levels obtained from Love et al (1077). 
Brailhwaite (1980). Murakami et al. (J999b), Dawson et al. (200(1). Pinter et al. (2002), Tindale and Roy (2002). and Van den heede et al (2(K)3). 
Alternative taxonomy, in brackets, taken from Ogle (1987). Distribution details wore compiled from various sources and. like ploidv level(s), will be 
dependent on taxon circumscriptions. Sample localities are only provided for New Zealand taxa, of which duplicate samples are identified by 
superscript letters. I lerbarium vouchers are given for samples from which novel sequence data was generated for this study. Sequences obtained from 
CenBank are referenced by the study that generated them. CenBank accession numbers are given for the tmL-tmF intergenic spacer and rhcL gene 
sequences included in the analyses. 


Taxon 

Ploidv 

«r 

Indigenous distribution: 
collection location for 
New Zealand samples 

Herbarium voucher/ 
previous study 

trnL-trnF 

intergenic 

spacer 

rbcL 
^ ;gene 

A . appendiculatum (Labi 11.) G.ProsI 

8x 

Australia. New Zealand: 

WELT P20492 

AY283202 


subsp. appendiculatum \A. termstrt * 


Mt. Cook 




Brownsey subsp. terrestre] 






A. appendiculatum subsp. 

8x 

New Zealand endemic: 

WELT P20493 

AY283203 


rnaritinmm (Brownsey) Brownsov 


Wellington 




| A. terrestre subsp. inaritimum 






Brownsey! 






A. bulbiftmun G.Forst subsp. 

4x 

New Zealand endemic: 

a WELT P20494 

aii AY283204 

a AY283226 

bulbiferum {A. bulbiferum 


A Palmerston North, 

"WELT P20495 



G.Forst., s.s.) 


B Stewart I. 




A. bulbiforum subsp. gracillimum 

8x 

Australia, New Zealand: 

a WELT P20496 

a AY283205 


(Coleuso) Brownsey (A. gracillimum 


A Hunterville, 

"WELT l , 20497 

"AY283206 


Colenso) 


H Stewart I. 




A. chnthanumse Brownsey 

4x 

New Zealand endemic: 

WELT P2049H 

AY283207 




Chatham 1. 




A . cimuwriorum Brownsey et He 

Hx 

New Zealand endemic: 

WELT P20499 

AY283208 


Lange 


Punakaiki 




A. flabellifoliuni (lav. 

4x, 5x, fix. 

Australia, New Zealand: 

WELT P20500 

AY283209 

AY283227 


8x, 12x 

Dannevirke 




A. flucciduin G.Forst subsp. flnc( idum 

4x 

Australia, New Zealand: 

a WELT P20501 

a,k AY283210 

a AY283228 

[A. flaccidum G.Forst s.s.) 


A Dunedin, H Paihia 

"WELT P20502 



A. flaccidum subsp. baurakiense 

4x 

N(!w Zealand endemic: 

WEI.T P20503 

AY283211 


Brownsov (A. haurakiense 


cultivated Auckland 




(Brownsey) Ogle) 







to 
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Tabu*; 1. Continued. 


Taxon 

A. hookcrianum Colenso 

A. lampmphyllum Carso 
A. lyullii (Hook.f.) I .Moore 
A. ohlonyifolium Colenso 

A. obtusntum (J.Korst suhsp. 
obtusutum [A. ohtustitum 
G.Forst., s.s.) 

A. ohtustitum suhsp. 

notihlandicum Hrownsey 
(/l. northlundicum 

(Hrownsey) Oglo) 

A. paupeivquitum 
Brownsov el l 1 .).Jacks. 

A. pulyodon G.Forst. 

A. richurdii (I look.f.) I look.I. 
A. scleroprium Homin'. 

A. shuttlt'uurthiunum kunzn 


Ploidv 


4 a 


4a 

Ha 

4a 


4a 


Indigenous distrihution: 
collection location for 
New Zealand samples 


Australia. New Zealand: x l)annevir 
"Banks Peninsula 
New Zealand endemic: Auckland 
New Zealand endemic: Castlepoint 
New Zealand endemic: A Paihia, 
"Palmerston North 
New Zealand, South America: 

A B1 uff. "Haast 


Ha 


Australia, New Zealand, 
Pacific Islands: Auckland 


Ha 


4 a 


Ha 

Ha 

Ha 


New Zealand endemic: 

Poor Knights Is. 

Asia. Australia, New Zealand, 
Pacific Islands: Palmerston North 
New Zealand endemic: Mt. Cook 
New Zealand endemic: Bluff 
New Zealand, Pacific Islands: 
cultivated Auckland 


A. trichomonas L suhsp. 
quadrivahms I l.K.Moy 
emend. Lovis 

A. tricliomanrs I.. suhsp. nov. 
A. acthiopicum (Burin.f.) Bech. 

A. antiquum Makino 
A . aurtnim (lav. 


(ex. Kermadnc Is.) 
4 a Sub-cosmopolitian: 

not available 



4a. Ha, Ha, 
10 a. 12a 


4 A 
4a 


Australia. New Zealand: Takaka 
Africa, Asia, Australia 

Asia 

Macaronesia 


x 


Herbarium voucher/ 
previous study 


A 


WELT P20504 
"WELT P20505 





*20506 
WELT P20507 
a WELT P20508 
"WELT P20509 
a WELT P20510 
"WELT P20511 


WELT P20512 


WELT P20518 

Vogel unpublished/ 
Pinter et al., 2002 
Yatabe et al. 2001 
Vogel unpublished 


trnL-trnF 

intergenic 

spacer 


rbcL 

gene 


A 


AY283212 
"AY283213 
AY283214 
AY283215 
A "AY283216 


AY283229 


AY283230 


AH 


AY28323 


A.H 


AY283217 


A "AY28323 


AY283218 


AY283224 

AF240666 


AY283236 

AF240654 


Al* 240657 


AB023502 

AF240642 


pc 




PC 


Z 


Z 




Z 


WELT P20513 

AY283219 

AY283233 


WELT P20514 

AY283220 

AY283234 

n 

WELT P20515 

AY283221 


s 

WELT P20516 

AY283222 


Sn 

WELT 1*20517 

AY283223 

AY283235 

* 

a . 




V m 

c 




Z 


to 4 


9 * % 









Table 1 . Continued 


Taxon 

A. austraiasicum (J.Sm) Hook. 
A. ceterach L. 

A. ensiforme Hook, et Grev. 

if 

A. hemionitis L. 

A. marinum L. 

A. nidus L. 

A. phillipsianum 
(Kiimmerle) Bir, 

Fraser-Jenk. 8: Lovis 


Ploidv 

Indigenous distribution: 
collection location for 
New Zealand samples 

4x 

Australia, New Guniea, 


Pacific Island 

2A', 4a, 6a 

Africa, Asia, Europe 

4a 

Asia 

? 

• 

Africa, Europe 

2a 

Africa, Europe 

4a 

Asia, Pacific Islands 

2a, 4a, 6a 

Africa 


A. prolongation Hook 

? 

• 

Asia 

A. sagittatum (DC.) Bange 

2x, 4 a 

Asia, Europe 

A. setoi N. Murak. et Seriz 

4a 

Asia 

A. theciferum (Kunth) Mett. 

4a 

Africa, America 

A. viride Huds. 

2a 

Africa, Asia, Europe, 



North America 

A. wrightii Hook. 

8a 

Asia 

1 iymenasplenium 

2a, ?4x 

Africa, Asia, Australia 

unilaterale (Lam.) Hayata 

%r 


Pacific Islands 



Herbarium voucher/ 
previous study 

trnL-trnF 

intergenic 

spacer 

rbcL 

gene 

WELT P20519 

AY283225 

AY283237 

Vogel unpublished 

AF240658 

AF240643 

Murakami et a/., 1999a 


AB014709 

Pinter et al ., 2002 

AF240663 

AF240848 

Pinter et al ., 2002 

AF240862 

AF240647 

'atabe et al. , 2001 


AB023500 

Listed as A. cord at um 

AF525235 

AF240650 


( rhumb.) Sw. by 
Pinter el al. (2002), but 
group with sequences 
identified as 
A. phillipsianum 
by Van den 
heede et aI. (2002). 


Murakami et al ., 1999a 


AR014891 

Vogel unpublished/ 

AF240661 

AF240848 

Pinter et al 2002 

Murakami et al ., 1999b 


AB013243 

Gastony and 


AF336099 

Johnson, 2001 

Pinter et al. , 2002 

AF240864 

AF240649 


Murakami et al ., 1999a 
Vogel unpublished/ 
Pinter et al 2002 


AF24066H 


AB014690 

AF240652 
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A. flaccidum 



Fig. 1 . Crossing polygon showing the Asplenium hybrid combinations known to occur naturally 
in New Zealand, updated from Brownsev (1977b, fig. 23). Thick line, intermediate line, dashed 
line, and no line, indicates hybridization is common, intermediate, rare, and unknown, re¬ 
spectively. Relative frequencies determined from analysis of 190 hybrid specimens cited bv 
Brownsev (1977b) and more recent collections in WELT. Subspecific taxa are not distinguished in 
the polygon because it is often difficult to determine which subspecies contributed to a hybrid. 


no doubt that hybridism plays an important part '. However, Brownsev (1977a, 
1977b) concluded that while hybridization amongst the species of New 
Zealand Asplenium was indeed relatively common, the spores of the hybrid 
plants were aborted and introgressant swarms were virtually absent. 

Brownsev (1977b) also noted that hybridization was unequally distributed, 
with some taxa participating frequently and others not at all. Although geo¬ 
graphy and ecology influence the frequency of hybridisation, some taxa are not 
known to form hybrids in New Zealand despite frequent svmpatric occurrence 
with congenerics. For instance, Asplenium polvodon , despite frequent svm- 
patrv, is not known to hybridize with A. bulbiferum subsp. bulbiferum , 
A. bulbiferum subsp. gracillimum , A. flaccidum subsp. flaccidum , or A. 
oblongifolium. Similarly, A. trichomanes, despite frequent sympatry. is not 
known to hybridize with A. bulbiferum subsp. gracillimum , A. hookerianum , 
or A. lyallii. 

A crossing diagram (Fig. 1) shows the extent and frequency of hybridization 
amongst New Zealand Asplenium taxa. Brownsev (1977b) suggested that the 
Asplenium taxa involved in hybridization in New Zealand were more closely 
related to one another, and constituted an 'Austral’ group of Asplenium. Spe¬ 
cies from this Austral group, comprising A. appendiculatum , A. bulbiferum , 
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A. chathamense, A. cimmeriorum, A. flaccidum, A. hookerianum , A. lamp- 
rophyllum, A. lyallii, A. oblongifolium, A. obtusatum, A. richardii, A. 
scleroprium and A. shuttleworthianum , are generally common in New Zea¬ 
land, and some of these species or their close relatives also occur in southern 
Australia, the southern Pacific, and South America. There are no diploid 
Asplenium species in New Zealand (Dawson et al., 2000), but amongst the 
tetraploid New Zealand representatives of the Austral Asplenium group. 
Brownsey (1977a, p.84) suggested the presence of “three natural groups of 
closely related species ... : (1) the A. flaccidum aggregate with thick and rather 
leathery bipinnate fronds, (2) the A. hookerianum!bulbiferum complex with 
thin highly dissected fronds, and (3) the A. lucidum!obtusatum complex with 
thick fleshy pinnate fronds" ( Asplenium lucidum G.Forst. non N.L.Burman is 
an earlier, but illegitimate, name for A. oblongifolium Colenso; Brownsey. 1979). 

Brownsey (1977b) considered the Asplenium species not hybridizing in New 
Zealand to have closer affinities with groups other than the Austral group. 
For instance, A. polyodon is widespread in the palaeo-tropics and has mor¬ 
phological similarities to species such as A. cuneatum and A. aethiopicum. 
Asplenium tnchomanes is widespread in the northern hemisphere. In Eu¬ 
rope, diploid, tetraploid and hexaploid forms of this species are known, and 
there is good cytological and morphological evidence that these hvbridize 
naturally with at least 13 other distinct species (Reichstein, 1981), in contrast 
to the total absence of hybridization in New Zealand. Asplenium flabellifolium 
also appears to have northern hemisphere affinities, showing some morpho¬ 
logical resemblance to A. viride. The relationship of A. pauperequitum is less 
clear. Brownsey and Jackson (1984) suggested that, based on morphology, this 
endemic species was not closely related to any of the other New Zealand 
species, although it may have distant affinities with A. polyodon. 

This mixture of postulated affinities suggests that the extant Asplenium taxa 
in New Zealand do not share a common New Zealand ancestor and that they 
are not monophyletic. Rather, it suggests they are derived from several non- 
New Zealand ancestors, implicating multiple, distinct origins of Asplenium in 
New Zealand. If multiple origins for New Zealand Asplenium are indeed the 
case, then ot considerable biogeographic interest is when these lineages ar¬ 
rived. Are the separate lineages ot New Zealand Asplenium old, perhaps in 
accord with the separation of New Zealand from the remainder of Gondwana 
about 80 million years ago (McLoughlin, 2001), and therefore consistent with 
a vicariant explanation? Or, do they appear younger than 80 m.y.a.. neces¬ 
sitating a long-distance dispersal explanation for their origins? 

Amongst the New Zealand members of the Austral Asplenium group there 
are eight tetraploid and nine octoploid taxa. Five of these octoploids are en¬ 
demic to New Zealand, and from morphological and ecological evidence they 

have been hypothesized to have arisen via auto- or allopolvploidy (Brownsev- 
1977b; Brownsey & de Lange. 1997). 

To conduct a phylogenetic analysis of the New Zealand Asplenium and to 
test vv hether their relationships are indeed reflected in hybridization frequency/ 
ability', we obtained DNA sequence data from two chloroplast markers. While 
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the chloroplast may only be inherited maternally, as has been reported by Vogel 
et al. (1D( in) tor European A. trichomanes, we anticipated that a chloroplast 
phylogeny of New Zealand Asplenium would nevertheless provide important 
insights into their origins, in terms of both the number of origins and, via 
a molecular clock approach, their timing. We also hoped to identify the 
chloroplast, or maternal, parent for each of the Austral octoploids. 

The trnL-trnF intergenic spacer was sequenced for every New Zealand 
Asplenium taxon presently recognized (Table 1), except A. trichomanes subsp. 
quaclrivalens for which we know of no extant population in New Zealand. 
We did not include in our study Pleurosorus rutifolius (R.Br.) Fee. another 
member of the Aspleniaceae indigenous to New Zealand, but its phylogenetic 
position is currently being investigated elsewhere (Johannes Vogel, pers. 
comm.). Not only does the trnL-trnF intergenic spacer amplify easily and 
consistently across a wide range of plants but, because it is non-coding, it 
evolves relatively fast, and accumulates both base-pair substitutions and 
insertion-deletion (‘indel) events that may be phylogenetically informative. It 
is therefore one of the best sequence markers available lor detecting genetic 
differences amongst the New Zealand Asplenium taxa. 

The rbcL gene was sequenced for representative New Zealand Asplenium 
taxa. Although the rbcL gene may evolve slower than the trnL-trnF intergenic 
spacer, it does not accumulate indel events, which means the alignment of se¬ 
quences from different taxa is easier. In addition, the large existing database of 
rbcL sequences for Asplenium taxa from around the world allows placement of 
representative New Zealand taxa into a global framework. 

Materials and Methods 

Table 1 details lor each of the indigenous Asplenium taxa currently rec¬ 
ognized from New Zealand (Brownsey and Smith-Dodsworth, 2000; Brownsey, 
2003) the taxonomy adopted for this study, the natural distribution, ploidy 
level, and the location of the sample(s) analyzed. We note that for some of 
these taxa there is debate as to whether they should be recognized at the 
subspecifc (Brownsey, 1977a) or specific (Ogle, 1987) level, and provide the 

alternative taxonomy where appropriate. Herbarium abbreviations follow 
Holmgren et al. (1990). 

DNA was extracted from silica-gel dried tissue using a modified CTAB pro¬ 
tocol (Doyle & Doyle, 1990). PCR was performed in 20 gL volumes containing 
lx Qsolution (Qiagen), lOmM Tris-HCl pH 8.8, 50 mM KC1, 1.5mM MgCl 2 . 250 
pmol dNTPs, 10 pmol of each primer, 1 U of Taq DNA polymerase (Qiagen). 
and approximately 50 ng of template DNA. 

The complete sequence for the trnL-trnF intergenic spacer was obtained 
Irom every Asplenium taxon indigenous to New Zealand (except Asplenium 
trichomanes subsp. quaclrivalens'. Table 1) using the primers E and F of 
Taberlet et al. (1991) and a thermocycling profile beginning with an initial de- 
naturation ot 94 C for 2 minutes, followed by 38 cycles of 94 C for 1 minute. 
58 C for 1 minute, and 72 C for 1 minute, with a final extension of 72°C for 
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5 minutes. Duplicate samples of the geographically-widespread tetraploid 
taxa in New Zealand [A. bulbiferum subsp. bulbiferum, A. hookerianum, A. 


flaccidum subsp. flaccidum, A. oblongifoli 


and A. obtusatum subsp. 


obtusatum), which are the most likely to be genetically variable, were analyzed 
to assess infra-taxon variation in this marker. Duplicate samples of the wide¬ 
spread octoploid A. bulbiferum subsp. grcicillimum were also analyzed. Local¬ 
ity details, herbarium voucher number, and GenBank accession numbers are 
given in Table 1. Sequence for the trnL-trnF intergenic spacer was also ob¬ 
tained from a cultivated plant of A. australasicum (J.Sm.) Hook. (Table 1). 

Sequences for the rbcL gene were obtained for representative New Zealand 
species from each of the three Austral groups, each of the non-Austral New 
Zealand species, and a cultivated plant of Asplenium australasicum (Table 1). 
The external primers aF and cR of Hasebe et al. (1994) were used with a 
thermocycling profile beginning with an initial denaturation of 94 C for 
2 minutes, followed by 38 cycles of 94°C for 1 minute, 58°C for 1 minute, and 
72°C for 1 minute, with a final extension of 72°C for 5 minutes. The novel 
internal primers rbcLAsForward2 (5'-AAGCCAAAATTAGGTCTATCTGC-3') 
and rbcLAsReverse2 (5-GGGAATTCTCTCGCAAAAAGAG-3) were used 
where necessary to obtain single amplification products and/or complete bi¬ 
directional sequencing. 

PCR products were purified using the CONCERT Rapid PCR Purification 
System (Gibco BRL). The purified PCR products were sequenced in both 
directions using an Applied Biosystems 373A DNA Sequencing System and 
the ABI PRISM™ Dve Terminate Cycle Sequencing Ready Reaction Kit 


(Perkin Elmer). 

Additional trnL-trnF and rbcL sequences were obtained from GenBank 
for species representative of the genetic diversity previously reported in the 
Aspleniaceae (Murakami et al., 1999a; Yatebe et al., 2001; Gastony and 
Johnson, 2001: Pinter et al., 2002; Van den heede et al., 2003). Details are 
provided in Table 1. Only rbcL sequence was available for some of the selected 
species, preventing their inclusion in the trnL-trnF alone and combined 
analyses (see below). 

In this study we follow Murakami (1995) in the recognition of Hymenasple- 
nium Hayata as a separate genus sister to Asplenium, and use H. unilateral, 
for which both trnL-trnF and rbcL sequences are available from GenBank, as 
the outgroup. However, to test the effect of outgroup choice (Adachi and 
Hasegawa, 1995), the rbcL analyses were repeated with H. unilaterale replaced 
with either H. hondoense (N. Murak. et Hatanaka) Nakaike (GenBank 
ABO14705) or H. riparium (Liebm.) N. Murak. (ABO14708). Previous studies 
(e.g., Murakami et al., 1999a; Van den heede et al., 2003) indicate these spe¬ 
cies are representative of the diversity known within Hymenasplenium. 

The sequences for each marker were aligned with ClustalX 1.8 (Thompson 
et al., 1997). The alignment of indels amongst the trnL-trnF sequences was 
edited further with Se-Al vl.O (Rambaut, 1995). One central region of the trnL- 
trnF intergenic spacer could not be unambiguously aligned between the {New 
Zealand Austral Asplenium species -I- A. australasicum } and the remaining 
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species. It was excluded from all analyses except those encompassing only the 
New Zealand Austral species and A. austmlasicum. Analyses of the trnL-trnF 
sequences were performed with sites encompassing indels either completely 
excluded, or included and treated as missing data. In addition, analyses were 
also performed with the large indels of A. marinum, A. aethiopicum, A. 
trichomanes, A. hemionitis, and the A. flaccidum group (the MATHF indels) 
included and treated as missing data, but with all other indels excluded. 

PAUP* 4.bl0 (Swoliord, 2002) was used to perform heuristic search analyses 
under maximum parsimony (MP; with tree bisection-reconnection branch- 
swapping, and 100 replicates of random sequence addition), and maximum 
likelihood (ML; with tree bisection-reconnection branch-swapping, and ran¬ 
dom sequence addition). Appropriate models of evolution for ML were se¬ 
lected using Modeltest v3.06 (Posada and Crandall, 1998) with the hierarchical 
likelihood-ratio test. Analyses were conducted on the trnL-trnF and rbcL data 
sets separately, and with them combined into a single data set after testing for 
incongruence using a partition homogeneity test as implemented in PAUP* 
4.bl0. Bootstrapping was performed to assess confidence in the groups iden¬ 
tified by maximum parsimony (n = 1000) or maximum likelihood (n = 100: 
with the parameters estimated by Modeltest fixed). 

The timing of divergences within Asplenium were estimated with the pro¬ 
gram r8s v.1.60 (Sanderson, 2002), which offers methods with the advantage 
of allowing rates of evolution to vary throughout the tree, and has been applied to 
other pteridophyte groups (e.g., Lycopodiaceae, Wikstrom and Kenrick. 2001; 
Equisetum, Des Marais et al., 2003). The tree used for the r8s analysis was that 
estimated by maximum likelihood. The r8s analysis was implemented under 
penalised likelihood, using the Powell algorithm with a log-scale pen¬ 
alty function, and a smoothing value of 125 as determined by cross-validation. 
Following the recommendation of Sanderson (2002), more distant outgroups 
[Polystichum vestitum (G. Forst.) C. Presl, GenBank AF208395; Thelypteris 
acuminata (Houtt.) Morton, D43919; Grammitis tenella Kaulf., AF468198) 
were used to estimate the branch length of Hymenasplenium unilaterale 
before being discarded prior to the r8s analysis. Gonfidence intervals (mean ± 
standard deviation) were calculated through a bootstrapping procedure with 100 
replicates (Sanderson and Doyle, 2001). The analysis was calibrated temporally 
by assigning an age of 140 m.y. (million years) to the most recent common 
ancestor of what appears to be the most distantly related extant components of 
the Aspleniaceae, Hymenasplenium and Asplenium. This date follows Skog 
(2001) who reported Aspleniaceae fossils from the early Cretaceous. 


Results 

Only two instances of trnL-trnF sequence variation within a taxon were 
found. The two samples of Asplenium bulbiferum subsp. gracillimum differ 
at position 50. The two samples of A. hookerianum also vary at position 50, 
as well as at two other positions invariant in the other taxa investigated. 
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Because of insertions and deletions, the length of the trnL-trnF sequences 
ranged between 238 and 398 base pairs. The trnL-trnF alignment, containing 
in dels, comprised 393 positions tor the Austral group plus Asplenium 
australosicuin, and 428 positions for all taxa with the central region that 
could not be unambiguously aligned excluded. Maximum parsimony analysis 
of all trnL-trnF sequences strongly supported recognition of an Austral clade, 
with A. austrolasicum recovered as its sister (not shown). As depicted in Fig. 2, 
MP analysis of the trnL-trnF sequences of the Austral taxa, together with A. 
austrolasicum as an outgroup, indicated three strongly supported sub-clades 
within the Austral clade, more or less corresponding to the groups suggested 
by Brownsey (1977a; see above). 

The Flaccidum clade, with 98% bootstrap support (BS), includes Asplenium 
appendiculatum suhsp. appendiculatum, A. appendiculatum subsp. mariti- 
mum, A. chathamense, A. flaccidum subsp. flaccidum, A. flaccidum subsp. 
haurakiense, A. lamprophyllum, and A. shuttleworthianum. The Bulbiferum 
clade, with 98% BS, includes A. bulbiferum subsp. bulbiferum, A. bulbiferum 
subsp. gracillimum, A. cimmeriorum, A. hookerianum, and A. richardii. The 
Obtusatum clade, with 100% BS, includes A. lyallii, A. oblongifolium, A. 
obtusatum subsp. northlandicum, A. obtusatum subsp. obtusatum, and A. 
scleroprium. The Bulbiferum and Flaccidum clades are supported as sister 
clades with 85% BS. The inclusion/exclusion of indel sites had minimal ef¬ 
fect on these principal results. 

Incorporation of representatives from each of the three Austral sub-clades 

in analyses ol rbcL sequences with a broader sample set, including species 

for which trnL-trnF sequences were not available, also supported recognition 

of the Austral clade. The rbcL sequences of A. oblongifolium and A. obtusatum 

were identical, so only the former was included in the analyses. The alignment 

ol the rbcL sequences, with no inference of indels necessary, comprised 1191 
base pairs. 

Using the model selected by Modeltest (TrN + I T G; base frequency = 
[0.2698, 0.2215, 0.2461, 0.2626]; rate matrix = [ 1 , 6.1972, 1 , 1 , 11.4499]; 
Invariable sites proportion = 0.5235; Gamma shape parameter = 0.8450), ML 
analysis selected a single tree with a —In likelihood score of 4426.3070. This is 
presented in Fig. 3, where bootstrap support is also indicated. Eight trees of 
487 steps (Cl = 0.624, RI = 0.689, RC = 0.430) are recovered under maximum 
parsimony. The consensus of these equally most parsimonious trees (not 
shown) is similar to the set of relationships depicted by the maximum like¬ 
lihood tree, except that the position of Asplenium antiquum within the 
Greater-nidus group is unresolved, A. bemionitis and A. phillipsianum form 
a clade that is sister to the remainder of the Marinum group, and A. wrightii is 
part of the basal polytomv rather than sister to the remainder of Asplenium. 

Pour major, monophyletic groups within this set ot Asplenium are evident, 
which for communication purposes we here label as the Aethiopicum, 
Ceterach-Phvllitis, Marinum, and Greater-nidus groups (Fig. 3). Asplenium 
wrightii does not appear to be closely related to any ol these groups, and it may 
represent a fifth group within Asplenium. Support for these four major groups 
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as A. theciferum and A. prolongatum. The Greater-nidus group is supported 
with 100% BS under MP and ML. The Austral group is supported with about 
80% BS under MP and ML. Within the Austral group, the Obtusatum (rep¬ 
resented by A. oblongifolium ) and Bulbiferum groups are supported as sister 
groups with about 70% BS under both MP and ML, which is in conflict with 
the trnL-trnF analysis. The species or group of species within the Greater-nidus 
group most closely related to the Austral group is not well resolved. 

Of the New Zealand species not belonging to the Austral group. Asplenium 
polyodon falls within the Aethiopicum group, and A. pauperequitum, A. 
flabellifolium, and A. trichomanes within the Marinum group. Asplenium 
flabellifolium and A. pauperequitum are strongly supported as sister species 
whereas A. trichomanes appears closest, of the species sampled, to A. viride. 
Use of Hymenasplenium hondoense or H. riparium as the outgroup instead of 
H. unilateral had little effect on these results, except in both cases BS for the 
Marinum group was less than 50% under MP. 

Gombination of the rbcL and trnL-trnF data sets is supported by a partition 
homogeneity test (P=0.14). Analyses of this combined data set are consistent 
with those of the rbcL data set in that four major, monophvletic groups can be 
recognized within the Asplenium sample set analysed here, and that the Austral 
group is part of the Greater-nidus group. With the MATHF indel sites included 
and treated as missing data but all other indel sites excluded, five most equally 
parsimonious trees of 659 steps (CI = 0.716, RI = 0.712, RC = 0.510) are recovered 
(not shown). These differ in whether the Obtusatum group or the Flaccidum 
group diverge most basallv within the Austral group, and whether the 
Aethiopicum group is sister to the Marinum group, an amalgamation of the 
Greater-nidus and Ceterach-Phvllitis groups, or to all three of these other major 
groups. Using the model selected by Modeltest (TIM + G; base frequency = 
[0.2777,0.2162, 0.2367,0.26941; rate matrix = [1,3.4608, 0.4752,0.4752. 5.12921; 
Gamma shape parameter = 0.3029), ML analysis selected a single tree with a—In 
likelihood score of 5766.1590 (not shown). This has the Flaccidum group 
diverging most basallv within the Austral group, and the Aethiopicum group 
sister to an amalgamation of the Greater-nidus and Ceterach-Phvllitis groups. In 
both instances the relationships are supported with less than 50% BS. 

A difference between analyses of the data sets is the sister group of the 
Greater-nidus group. In the rbcL data it is the Aethiopicum group (61% BS MP. 
76% BS ML), w^hereas in the combined rbcL and trnL-trnF data set it is the 
Ceterach-Phvllitis group (93% BS MP, 97% BS ML; the Greater-nidus and 
Ceterach-Phvllitis sister relationship is also recovered in analysis of the trnL- 
trnF data set alone). In addition, the bootstrap support for the Marinum group 
is low^er in the combined rbcL and trnL-trnF data set (<50% BS MP, 56% BS 
ML) than the rbcL data set (70% BS MP, 79% BS ML). When all indel sites in 
the combined rbcL and trnL-trnF data set are excluded the four major groups 
are still recovered (not shown), albeit with the Marinum group with onlv low 
BS, but their relationships to each other are unresolved. 

The age of the most recent common ancestor of the Greater-nidus group w r as 
calculated by r8s, under penalised likelihood with the Pow r ell algorithm and 
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a log-scale penalty function, at 34 ± 5 million years ago (m.y.a.). Similarly, the 
ages calculated by r8s for the most recent common ancestor between each of 
the non-Austral New Zealand species and their closest non-New Zealand 
relative in this data set were: Asplenium polyodon and A. ensiforme, 31 ± 6 
m.y.a.; {A. flabellifolium + A. pauperequitum} and A. mcirinum, 43 ± 7 m.y.a.; 
and A. trichomcines and A. viride, 30 ±7 m.y.a. Substituting H. hondoense or 
H. riparium as the outgroup instead of H. unilaterale had little effect on these 
results, producing age estimates differing by only ±2 m.y. 

The trnL-trnF sequences of the tetraploids Asplenium flaccidum subsp. 
flaccidum and A. chathamense . and the octoploid A. appendiculatum subsp. 
maritimum are identical to one another. The trnL-trnF sequences of the 
tetraploid A. flaccidum subsp. haurakiense and the octoploid A. appendicu¬ 
latum subsp. appendiculatum are also identical, and differ from the three 
aforementioned taxa by a single shared base-pair substitution, at position 50. 
The trnL-trnF sequence of the octoploid A. shuttleworthianum differs from all 
others sampled by several base-pair substitutions and indel events. 

The octoploids Asplenium bulbiferum subsp. gracillimum and A. cimmer- 
iorum share several putative synapomorphies in their trnL-trnF sequence with 
the tetraploid A. hookerianum . and differ from the tetraploid A. bulbiferum 
subsp. bulbiferum by between five and seven base-pair substitutions. Depend¬ 
ing on the reconstruction of character evolution, the octoploid A. richardii 
shares two or three trnL-trnF sequence putative synapomorphies with A. 
hookerianum, A. bulbiferum subsp. gracillimum, and A. cimmeriorum, but 
branches basally to them and is sister to that clade of species. 

The trnL-trnF sequences of the octoploids A. lyallii and A. scleroprium are 
identical to those of the tetraploids A. oblongifolium and A. obtusatum, while 

that of the octoploid A. obtusatum subsp. northlandicum exhibits two 
autapomorphies. 


Discussion 

The evolutionary relationships inferred here from chloroplast sequences 

support the hypothesis based on hybridization frequency that New Zealand 

Asplenium are non-monophyletic, and that there is one large, relatively closely 

related ‘Austral’ group to which several other New Zealand species are only 

distantly related. Members of the Austral group hybridize with one another in 

New Zealand, but not with the non-Austral species. Similarly, Murakami et al. 

(1999a) observed natural hybridization to only occur between closelv related 
Asplenium species. 

Our analyses indicate at least four major monophyletic groups within the set 
of Asplenium investigated here, and our findings more or less correspond to 
those of other studies of chloroplast sequences (eg. Murakami et aL, 1999a; 
Gastonv and Johnson, 2001; Pinter et al., 2002). However, while the support in 
our analyses for each of these groups is generally strong (the Marinum group 
being an exception in some analyses), the relationships between them are not 
clear. This may be resolved by additional sampling, which may also indicate 
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the presence ot other major monophyletic groups (e.g., A. wrightii?) or a re- 
circumscription of the groups outlined here. 

Species from the Austral group have not previously been sampled genet¬ 
ically. Although Schulze et al. (2001) reported a rbcL sequence (GenBank 
AF318601) for a plant of 'A. bulbiferum'’ from the Heidelberg Botanical 
Garden, their sequence is very different from that of the wild New Zealand 
material of A. bulbiferum reported here. Of the species here sampled. 
AF318601 appears closest to A. theciferum (not shown). The placement of 
the Austral group as a well-supported sub-group within the Greater-nidus 
group is perhaps unsuspected horn morphological comparison. 

The Greater-nidus group itseli appears to have a predominantly southern 
hemisphere and north-western Pacific distribution, with Asplenium prolonga- 
tum also extending to India and China. Available chromosome counts (Love 
et al., 1977; Murakami et al. , 1999b; Dawson et al., 2000; Tindale and Roy. 
2002) indicate that all of the species sampled here from the Greater-nidus 
group are at least tetraploid, suggesting that this may be the ancestral state of 
the group. Further sampling and investigation is required to confirm this. Also 
ot interest in the Greater-nidus group is the apparent lack of a close 
relationship between A. antiquum and the other birds' nest ferns, or between 
the Austral group and the other members of the Greater-nidus group with 
divided laminae (i.e., A. prolongatum and A. theciferum). 

Even with this limited sampling, it is clear that New Zealand’s Asplenium 
flora has had multiple independent origins from distantly related parts of the 
genus. There has been at least one separate migration in the Greater-nidus 
group, at least one in the Aethiopicum group, and at least two in the Marinum 
group. Moreover, molecular dating clearly indicates that these origins can be 
attributed to dispersal, and not vicariance. All age estimates for the most recent 
common ancestor of each New Zealand species and their closest non-New 
Zealand relative in this sample set are younger than 80 m.y.a., which is 
approximately when the New Zealand landmass separated from Australia and 
the rest of Gondwana. The closest age estimate to this boundary is the 43 m.y.a. 
between {A. flabellifolium + A. pauperequitum) and A. marinum. However, 
this is probably, at least in part, reflecting a distant relationship between A. 
marinum and A. flabellifolium or A. pauperequitum. In all cases, inclusion 
of non-New Zealand samples more closely related to each New Zealand 
Asplenium would make the age estimates younger. Also consistent with 
a dispersal interpretation for the origins of New Zealand Asplenium. 
Mildenhall (1980) gives the mid-Miocene (about 15 m.y.a ) as the earliest ap¬ 
pearance of Asplenium spores in the New Zealand fossil record. 

Of considerable importance to the molecular dating is the use of an age of 
140 m.y. for the most recent common ancestor of Asplenium and Hvmenas- 
plenium to calibrate the genetic divergences within Asplenium. If the sep¬ 
aration between Asplenium and Hymenasplenium was actually older than 140 
m.y., then the divergences estimated here within Asplenium would also be 
older. However, constraining the most recent common ancestor of the Greater- 
nidus group to 80 m.v.a., which is the youngest age consistent with a vicariant 
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origin of the New Zealand Austral Asplenium species amongst this data set, 
results in a calculation of 325 m.y.a. for the divergence between Asplenium 
and Hymenasplenium. This is quite inconsistent with the fern fossil record 
(Skog, 2001), indicating that even allowing for errors in the calibration, 80 m.y. 
old vicariant origins for New Zealand Asplenium can be rejected. 

Alternatively, although the earliest known Aspleniaceae fossils may be ap¬ 
proximately 140 m.y. old, the separation between Asplenium and Hymenas¬ 
plenium could be more recent (and perhaps much more). If so, the divergences 
within Asplenium would also be younger than estimated here, but would still 
necessitate inference of dispersal rather than vicariance for the origins of 
New Zealand Asplenium. Indeed, an important caveat is that because the 
lossil record places, as tar as we know, no ‘younger’ bound on the divergences 
investigated here, their real ages could be orders of magnitude younger 
than calculated. 

Within the Austral group, the following taxa also occur outside New Zea¬ 
land: Asplenium appendiculatum subsp. appendiculatum, A. bulbiferum 
subsp. gracillimum, A. flaccidum subsp. flaccidum, A. hookerianum, A. 
obtusatum subsp. obtusatum, A. obtusatum subsp. northlandicum, and A. 
shuttleworthianum (Table 1). From the molecular dating analysis above, the 
most recent common ancestor ol all ot these taxa is younger than the separa¬ 
tion ot New Zealand from Australia. Consequently, inference of at least one 
dispersal event for each of these taxa is required to explain their occurrence in 
New Zealand and elsewhere. Even leaving aside the molecular dating evidence 
above, it seems unlikely, as pointed out by Brownsey (2001), that the disjunct 
populations of each ol these taxa have remained sufficiently unchanged mor¬ 
phologically to be regarded as the same taxon through some 80 m.y. of sep¬ 
aration, if thev were indeed vicariant. 

Recent molecular studies have found long-distance dispersal to be prevalent 
in the origins of the New Zealand flora (reviewed by Winkworth et al., 2002), 
with few exceptions (Stockier et al., 2002). Within ferns, molecular dating has 
suggested an origin via dispersal for New Zealand Polystichum (Perrie et al., 
2003a). In addition to receiving elements via long-distance dispersal. New 
Zealand has also acted as a source for other regions (Wright et al., 2000; 
Lockhart et al., 2001), and Brownsey (2001) suggested Asplenium bulbiferum 
subsp. gracillimum may have dispersed from New Zealand to Australia. 

The Austral group ot Asplenium is relatively widespread in the southern 
Pacific and Australasian regions. Conspecifics or close relatives of the New 
Zealand taxa from the Bulbiferum group are known to occur in Australia: the 
Flaccidum group in Australia and some south Pacific islands; and the 
Obtusatum group in Australia, South America and some south Pacific islands. 
Further sampling may indicate that species from outside the Australasian and 
southern Pacific regions, or other species present in these regions but not in 
New Zealand, may also belong to this Austral group. 

Although each appears strongly circumscribed, the relationships among the 
Bulbiferum, Flaccidum, and Obtusatum groups are not well resolved because 
of one of the few points of conflict between the two molecular data sets studied 
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here. The trnL-trnF data strongly suggests that the Bulbiferum and Flaccidum 
groups are sister groups, whereas almost equally strongly the rbcL data sug¬ 
gests that the Bulbiferum and Obtusatum groups are each others’ closest rela¬ 
tives. The underlying explanation for this conflict is not known, and it is likely 
that more character data will be required to resolve this issue. However, it can 
be noted that the highest frequency of inter-group hybridization is between the 
Bulbiferum and Flaccidum groups, with Asplenium bulbiferum subsp. bulbife¬ 
rum X A. flaccidum subsp. flaccidum and A. bulbiferum subsp. gracillimum X 
A. flaccidum subsp. flaccidum being particularly common (Fig. 1). 

Asplenium lamprophvllum, despite its creeping rhizome and broad lamina 
segments, falls within the Flaccidum group whose other constituents have 
erect rhizomes and narrow lamina segments. Also in the Flaccidum group is A. 
shuttleworthianum, which with its narrow lamina segments and marginal sori. 
bears a strong morphological resemblance to, and indeed can be difficult to 
separate from, the Pacific A. gibberosum (G.Forst.) Mett. The latter is some¬ 
times placed (e.g., Brownlie, 1977) within the segregate genus Loxoscaphe 
T.Moore, as L. gibberosum (G.Forst.) T.Moore. However. rbcL sequence 
(Gastony and Johnson, 2001) from African material of the type species of 
Loxoscaphe, L. thecifera (Kunth) T.Moore (= A. theciferum), does not fall 
within the Austral group (Fig. 3). If A. gibberosum , for which sequence data is 
presently unavailable, is more closely related to A. shuttleworthianum than 
A. theciferum, then at least some authors' (e.g., Brownlie, 1977) circumscrip¬ 
tions of Loxoscaphe are polyphyletic. Alternatively, if A. gibberosum is more 
closely related to A. theciferum than A. shuttleworthianum, such that 
Loxoscaphe sensu Brownlie (1977) is monophyletic (albeit nested within 
Asplenium : Gastony and Johnson, 2001), then the resemblance of A. shuttlewor¬ 
thianum to A. gibberosum represents another case of striking morphological con¬ 
vergence within Asplenium (Murakami et al., 1999a; Van den heede et al., 2003). 

Of the species investigated, the closest relative of Asplenium pauperequitum 
was found to be A. flabellifolium, and not A. polyodon as tentatively suggested 
by Brownsev and Jackson (1984). Asplenium pauperequitum is one of the 
rarest ferns in New Zealand, being known only from the Poor Knights Islands 
(Brownsey and Jackson, 1984). Nevertheless, the relationship between A. 
flabellifolium and A. pauperequitum is not close. These two species show little 
morphological similarity, and r8s, using the calibration described above, es¬ 
timates a divergence time between them of 19 m.y. (with the caveats indicated 
above). Further sampling, particularly in the Pacific region, may find species 
with greater affinity to A. pauperequitum, and possibly even its tetraploid 
progenitor(s). 

Brownsey (1977b) suggested that the octoploid Asplenium appendiculatum 
was probably derived from the tetraploid A. flaccidum via autopolyploidy. The 
discovery that the trnL-trnF sequence of A. appendiculatum subsp. appendi¬ 
culatum is identical to that of A. flaccidum subsp. haurakiense, while the trnL- 
trnF sequence of A. appendiculatum subsp. maritimum is identical to those 
of A. chathamense and A. flaccidum subsp. flaccidum may indicate two 
independent polyploid events, mirroring other instances where different chlo- 
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roplast sequences have been found in both the putative progenitors and 

descendants (Soltis and Soltis, 1999). While certainly deserving of further 

investigation, the trnL-trnF sequences of A. appendiculatum subsp. appendi- 

culatum and A. flaccidum subsp. haurakiense differ from the others only in 

having a thymine rather than an adenine at position 50. This base pair appears 

particularly prone to homoplasy, with reversion between an adenine and 

thymine occurring independently at least once in the Flaccidum group, at least 

twice in the Bulbiferum group, and again at least once outside the Austral 

group. The propensity tor change shown by position 50 mav be due to its 

central location in what appears, if transcribed, to be a hairpin loop-forming 

stretch of the trnL-trnF intergenic sequence. An indication of two independent 

polyploid origins for A. appendiculatum from these data is dependent on the 

associated change at position 50 not being homoplastic, which it clearly is in 
other instances. 

The chloroplast sequence of Asplenium shuttleworthianum is quite distinct 

from any of the other New Zealand taxa sampled here from the A. flaccidum 

group, (xiven the primarily tropical distribution ol A . shuttleworthianum , 

its maternal tetraploid progenitor, il still extant, is probably to be found in 
the Pacific. 

Intriguingly, the trnL-trnF sequences from the octoploid Asplenium bulbife- 
ruin subsp. gracillimum implicate A. hookenanum as its chloroplast or maternal 
parent rather than A. bulbiferum subsp. bulbiferum . This raises the possibility 
that A. bulbiferum subsp. gracillimum is an allopolvploid between the non- 
bulbiferous A. hookenanum , from which it has inherited its chloroplast, and 
A• bulbiferum subsp. bulbiferum , with which it shares a close morphological 
similarity including bulbil production. This finding might be taken as support 
for the recognition at species lev el ot the two entities here regarded as subspecies 
of A. bulbiferum (Ogle, 1987; Table 1), and this is presently being investigated 

in prep.). The polymorphism shared between 
A. bulbiferum subsp. gracillimum and A . hookerianum is at position 50 which. 

because ot the suspicion outlined above that this base pair is susceptible to 

homoplasy, should be interpreted cautiously in any inference of multiple 
polyploid events. 

Asplenium cimmeriorum has an identical trnL-trnF sequence to one of the 

A * bulbiferum subsp. gracillimum samples. Whether the non-bulbiferous A. 

cimmeriorum is similarly derived from an allopolyploid event between 

A. bulbiferum subsp. bulbiferum and A. hookerianum — either directly via an 

independent event or following divergence from A. bulbiferum subsp. 

gracillimum — or is an autopolyploid ot A. hookerianum requires further 

research of the sort employed by Perrie et al. (2003b) using genetic markers 

from throughout the genome. Asplenium richardii is sister to a clade that 

includes A. hookerianum , from which Brownsey (1977a) suggested the former 
was possibly an autopolyploid. 

Brownsey (1977b) regarded Asplenium lvallii and A. scleroprium as prob¬ 
able allopolyploids based on their close morphological resemblance to sterile 
hybrids between members of the Bulbiferum and Obtusatum groups, and 
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Flaccidum and Obtusatum groups, respectively. The trnL-tniF sequences 
clearly implicate a chloroplast or maternal parent from the Obtusatum group 
for both A. lycillii and A. scleroprium. However, which of the two extant 
tetraploid species was involved in either case cannot be ascertained. Despite 
considerable differences in their spore morphology (Brownsey, 1977a), A. 
oblongifolium and A. obtusatum share identical sequences for both the trnL- 
trnF and rbcL markers. Their psbC-trnS sequences, another chloroplast 
intergenic spacer, are also identical (unpub. data). In contrast, the trnL-trnF 
sequence of A. obtusatum subsp. northlandicum exhibits two autapomorphies 
relative to A. obtusatum subsp. obtusatum , from which Brownsey (1977b) 
suggested the former arose by autopolyploidy. 

New Zealand appears to have been colonized, via dispersal, by several major 
groups within Asplenium. Most of the New Zealand species, however, belong 
to the Austral sub-group of the Greater-nidus group, and their close relation- 

to one another is reflected in their ability to hybridize. Future work will 
involve obtaining sequence data from Australian and Pacific material to assess 
the relationships and migration patterns within Asplenium in this region. 
Investigation into the relationships amongst the Asplenium species in New 
Zealand will also continue, with an emphasis on understanding the evolu¬ 
tionary history of some of the polyploids. 
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